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The Neutron Electric Dipole Moment in the Instanton Vacuum: 
Quenched Versus Unquenched Simulations 
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We investigate the role played by the fermionic determinant in the evaluation of the CP-violating 
neutron electric dipole moment (EDM) adopting the Instanton Liquid Model. Significant differences 
between quenched and unquenched calculations are found. In the case of unquenched simulations 
the neutron EDM decreases linearly with the quark mass and is expected to vanish in the chiral 
limit. On the contrary, within the quenched approximation, the neutron EDM increases as the 
quark mass decreases and is expected to diverge as 1/m f in the chiral limit. We argue that such 
a qualitatively different behavior is a parameter-free, semi-classical prediction and occurs because 
the neutron EDM is sensitive to the topological structure of the vacuum. The present analysis 
suggests that quenched and unquenched lattice QCD simulations of the neutron EDM as well as 
of other observables governed by topology might show up important differences in the quark mass 
dependence, for m q < Aqcd- 



1. INTRODUCTION 

The electric dipole moment (EDM) of the neutron pro- 
vides direct information on the violation of the parity and 
time-reversal symmetries. Both the strong and the elec- 
troweak sectors of the Standard Model (SM) can gener- 
ate violations of the above symmetries. As for the strong 
sector, it is known [1| that a gauge-invariant definition 
of the QCD vacuum requires to supplement the classical 
action with an additional gauge-invariant and renormal- 
izable term, which in Euclidean space reads 
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is the topological charge operator, F^ v — (1/2) e^uapFap 
is the dual gluon field strength (incorporating the strong 
coupling constant) and 9 is a (real) dimensionless param- 
eter. The term (11.2(1 is a source of CP violation and goes 
under the name of strong 0-term. 

A second independent physical origin for a term of the 
form Q1.2[l comes from the weak sector of the SM. The 
observation of CP violation in K-meson systems implies 
that the quark mass matrix M is not real and the mass 
term in the Lagrangian has the general form 



The mass matrix can be made real and diagonal by means 
of an appropriate chiral rotation, which generates a shift 
in the 9 parameter: 



9 



9 + argdet(M) 



(1.5) 



C 



-rR 



M 



^ Ma V L + & M/, ^ R 



(1.4) 



The real constant 9 is an additional dimensionless pa- 
rameter, which has to be fixed from experiment. This can 
be done by exploiting the fact that (|1.2|l is a source of 
CP violation which leads to a non- vanishing value of the 
neutron EDM. At present there are only upper bounds 
on the neutron EDM and the most constraining one is 
|d„| <6.3x 1CT 13 (e-fm) @|. 

In order to translate this experimental information 
into a constraint on 9, one needs to compute the neu- 
tron EDM in QCD, including the contribution of the 
topological term i9 Q. So far, this has been done only 
within model-dependent frameworks, starting from the 
works of Refs. y, 3 In Q Baluni performed a calcu- 
lation of the EDM in the Bag Model and found |d n | = 
2.7 x 1CT 3 9 (e • fm). In QJ Crewther et al. proposed an 
approach based on current algebra relations and found 
|d„| =3.6x 10" 3 6 (e-fm). 

The above model calculations agree in pointing out 
that \9\ < 2 • 10~ 10 . Understanding why 9 is so small is 
an open challenge, which goes under the name of strong 
CP problem. 

In order to estimate the relevant matrix element in 
a model-independent way, a non-perturbative approach 
based on the fundamental theory, like lattice QCD, is re- 
quired. However a lattice estimate of the neutron EDM 
is not yet available Q. Recently [6j a new strategy for 
computing the neutron EDM in lattice QCD has been 
proposed. The starting point is the expansion of the 



matrix element of the EDM operator to lowest order in 
9. This allows to remove the complex #-term (|1.2|) from 
the action, but at the expense of calculating the matrix 
element of a composite operator made of the electromag- 
netic dipolc and the topological charge operators: 
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where \N(8)) denotes the neutron state at a finite value 
of 9. 

Computing the matrix element (J1.7JI on the lattice is 
a very challenging task. The main source of difficulties 
is that the topological charge operator is very noisy . 
A possible way-out, suggested in Refs. y| 0- consists 
in exploiting the anomalous Ward identities to replace 
the topological charge operator in l|1.7|) with the pseudo- 
scalar density operator, namely 
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where ( ) denotes the quantum average over all config- 
urations, Nf is the number of quark flavors, mo is the 
bare quark mass (1/mo = (l/m u + \/m c i + l/m s )/3 for 
Nf = 3) and P is the flavor-singlet pseudo-scalar opera- 
tor P(x) = T,f^f( x h5ipf{x). In the r.h.s. of Eq. (fT£|) 
( )disc denotes the quantum average obtained by retain- 
ing only the Wick contractions in which the pseudo-scalar 
operator is contracted in a virtual quark loop. The nu- 
merical applicability of such a strategy is still under in- 
vestigation. 

The use of the quenched approximation is quite nat- 
ural for a first-time calculation. We point out however 
that neglecting the fermionic determinant may result in 
a sizable systematic error in case of the neutron EDM, 
since the evaluation of the latter involves flavor-singlet 
operators. Indeed, in case of hadron masses and decay 
constants quenched calculations turn out to be accurate, 
since they appear to agree with available (partially) un- 
quenched simulations within 5 -r- 10% accuracy ja |j|. 
However, in case of observables which are directly related 
to the topological properties of the vacuum, the contri- 
bution of the fermionic determinant may become more 
important. 

This can be seen, for example, by considering the Index 
Theorem, according to which the total topological charge 
of a gauge configuration relates to the difference of the 
number of left- and right-handed zero-modes of the Dirac 
operator: 



On the other hand, the fermionic determinant can be 
written in terms of eigenvalues of the Dirac operator as: 
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where Ipipx = A^X) an d v is the number of zero-modes. 

Equations l|1.9fl and (| 1 . 1 U|) imply that, for small val- 
ues of the bare quark masses, the fermionic determinant 
suppresses configurations with a non- vanishing topologi- 
cal charge. Since the EDM is an observable which relates 
directly to topology, we argue that evaluating this quan- 
tity in quenched and unquenched lattice QCD might lead 
to quite different results. Of course, one can still hope 
that sufficiently away from the chiral limit quenched and 
full simulations give comparable results. 

While waiting for lattice results on the neutron EDM, 
the aim of this work is to investigate the role played 
by the fermionic determinant using a model which is 
expected to take into account properly topological ef- 
fects. Such a model is the Instanton Liquid Model 
(ILM), since instantons are topological gauge configu- 
rations which dominate the QCD Path Integral in the 
semi-classical limit. They generate the so-called 't Hooft 
interaction, that solves the U(l) problem jj and sponta- 
neously breaks chiral symmetry [l(j, but does not con- 
fine. Evidence for such an instanton-induced interaction 
in QCD comes from a number of phenomenological stud- 
ies [ll| , as well as from lattice simulations [l^, HjJ U_ll • 
The ILM assumes that the QCD vacuum is saturated 
by an ensemble of instantons and anti-instantons. The 
only phenomenological parameters in the model are the 
instanton average size and density. Their values were es- 
timated long ago from the global vacuum properties [l5| . 
Since this model provides very successful descriptions of 
both the pion and nucleon electromagnetic structure |l6( 
and of the topological properties of the QCD vacuum 
fl7| , we expect its prediction for the neutron EDM to be 
realistic. 

The dynamical mechanism leading to the formation of 
the neutron EDM in the instanton vacuum was recently 
investigated by one of the author 18]. It was found that, 
during the tunneling processes, the 0-term generates an 
effective repulsion between matter and anti-matter, in 
the neutron. As a consequence, quarks and anti-quarks 
migrate in opposite directions. Hence, at least on the 
semi-classical level, the EDM arises from the local sep- 
aration of positive and negative baryonic charges in the 
neutron. It does not follow from the displacement of the 
positive and negative electric charge carried by the va- 
lence quarks, as one would intuitively expect in a naive 
non-relativistic quark model picture. 

In this work we show that the use of the ILM clearly 
suggests that neglecting the fermionic determinant leads 
to a dangerous divergence in the chiral limit, which in the 
full model is regulated by topological screening. We will 
compare unquenched and quenched model calculations at 
values of the space-time volume and of the current quark 
mass which are comparable to those used in present-day 



lattice simulations, and provide numerical estimates of 
the neutron EDM. 



2. RELATING THE NEUTRON EDM TO 
EUCLIDEAN CORRELATION FUNCTIONS 



In order to model quenched and full QCD, we will use 
the Interacting Instanton Liquid Model (IILM) developed 
by Shuryak and collaborators (see Ref. [lj for a review) . 
In this model the configurations of the instanton ensem- 
ble are generated by means of a Metropolis algorithm 
which explicitly accounts for the fermionic determinant. 
The instanton size distribution is not fixed a priori, but 
it is obtained dynamically. The average density was fixed 
by minimizing the free energy [19J. In the IILM the ne- 
glect of the fermionic determinant generates pathologies 
which are analogous to the ones observed in quenched 
lattice QCD (see Ref. [l4|). 

The feasibility of instanton calculations relies on two 
important simplifications, which arise when one restricts 
the functional integration over the gauge field configu- 
rations to an ordinary integral over the collective coor- 
dinates of an ensemble of iV+ instanton and N— anti- 
instantons. On the one hand, one has to deal with a dra- 
matically smaller number of degrees of freedom, typically 
of the order of 12 x (JV + + iV_) ~ 10 3 (for a typical box of 
the size of those used in lattice simulations). On the other 
hand, evaluating the topological charge operator in the 
instanton vacuum is trivial, because it amounts to simply 
counting the number of instantons and anti-instantons 
in the ensemble (indeed Q = N+ — N- ). Due to these 
simplifications, the numerical simulations leading to the 
EDM can be performed on a regular workstation. 

We show that, within the IILM, quenched simula- 
tions become unreliable for quark masses smaller than 
the strange quark mass. Indeed, for m q > 200 MeV, 
quenched and full simulations give results which agree 
within statistical errors, whereas for quark masses of the 
order of 100 MeV, quenched calculations overestimate 
the EDM by a factor ~ 4. This is due to the appearance 
of the chiral divergence mentioned previously. Moreover, 
we stress that other quantities, like e.g. the nucleon mass 
or the quark condensate, are not drastically affected by 
quenching. 

Assuming a linear mass dependence, the extrapolation 
of the unquenched ILM results to a physical value of the 
light quark mass between 4 and 10 MeV yields: |d„| = 
(6 4- 14) x 10~ 3 6 (e • fm), which is a factor about 2-^4 
larger than the estimates obtained in Refs. [3,l3l- 

The paper is organized as follows. In Section[3we shall 
discuss our method for relating the neutron EDM to Eu- 
clidean correlation functions. In Section we will study 
these correlators in the instanton vacuum and derive the 
quenched and unquenched expressions for the neutron 
EDM. In Section 0] we analyze the chiral behavior both 
for quenched and unquenched calculations. Results are 
presented and discussed in Section [5J and, finally, con- 
clusions are summarized in Section [5] 



To compute the neutron EDM in a field-theory frame- 
work, we start by considering the following Euclidean 
correlation function: 



Z^'(p,q) = |rfx|dye ipx+iqy 2/3G7'(x,y) 
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dq 3 
G^'(x,y) = (0|T[Jft(0 J 2r)J 4 (y,T)^(x,0)|0) (2.2) 

where Jji(x) = e a bc{u^C^db)dJ is an interpolating op- 
erator for the neutron, p = (p, oj), p' = (p + q, U)'). In 
the limit of large Euclidean times r one can show that: 
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where M is the mass of the neutron, Z is a constant 
defined by 



(0\J N (0)\N(p, S )) = Zu(p, S ), 
where u(p, s) is a Dirac spinor, and 
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with 043 the Euclidean version of O03 = *7o73- Notice 
that the form factors Fa(q 2 ) and F 3 (q 2 ) do not appear 
if we assume that the matrix element is invariant under 
C, P, and T. In particular F 3 (q 2 — 0) relates directly to 
the EDM through the relationship: 
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In order to isolate the contribution to the EDM it is 
convenient to define the following trace: 



E T (p,q) = Tr[£ r (p,q) T 3 ], E3 = ( ^ ^ 

and obtain: 
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Taking q = we obtain: 
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numerical integrations can be further reduced to two, 
by exploiting the rotational symmetry of G r (x = 0,y) 
around the y 3 axis: 

£ r (x = 0,q = 0) = 2tt f dy 3 y 3 [ d\y±\ \y±\ 

J —OO J 

■ Tr[G T (x = 0,y)r 3 ] . (2.13) 

As usual, the neutron mass M and the coupling Z 
to the neutron interpolating operator appearing in the 
r.h.s. of Eq. Ij2.12|) can be extracted from the zero mo- 
mentum two-point correlator, 
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Fig. 1: ILM results for log(Tr[P T (p = 0) 74]). The linear be- 
havior denotes the isolation of the nucleon pole. The slope of 
the line corresponds to the nucleon mass M (in units of inverse 
fm) while its intercept corresponds to log (2 Z ). These results 
were obtained in the quenched approximation, by averaging 
over the configurations of an ensemble of 56 pseudo-particles 
in a volume V = 2.4 3 x 4 = 56 fm 4 and quark masses of 
m q — 1.00 A.qcd- 



All correlation functions defined so far are projected 
onto fixed momentum components, hence the spatial po- 
sition of the sources and of the electromagnetic operator 
are completely unspecified. From the practical point of 
view, evaluating these quantities can be very challenging, 
because they involve a six-dimensional integration over 
the two spatial hyperplanes. This has to be performed 
numerically and requires to compute a large number of 
evaluation points for the integrand, each of which is a 
three-point correlation function. However, as long as we 
are interested in a static property of the neutron, such 
as the EDM, it is possible to reduce the number of spa- 
tial integrations by working in a mixed space-momentum 
representation, where one projects only on the momen- 
tum transfer by the photon and leaves the source and 
the sink at a fixed spatial location. In this way one gives 
up the information about the initial momentum of the 
neutron. Let us define the correlator: 
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Choosing x = we obtain: 
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Evaluating S T (x = 0, q = 0) involves only a three- 
dimensional integration over the position of the electro- 
magnetic current. Moreover, the number of independent 



3. THE NEUTRON EDM IN THE INSTANTON 
VACUUM 

In QCD the neutron EDM is non-zero only at finite 
values of the 0-angle. To lowest order in this parameter 
the three-point correlation function l|2.2|l reads: 
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where |0) refers to the vacuum at 9 = 0. 

In the instanton vacuum the topological charge is con- 
densed around instantons and anti-instantons, and corre- 
lation functions are computed by averaging over the con- 
figurations of a grand-canonical ensemble of such pseudo- 
particles with a partition function given by: 
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In this formula Si n t is the bosonic instanton-instanton 
interaction, z% denotes the position of the pseudo-particle 
of size pi, dUi is the Haar measure normalized to unity 
and d(pi) is the instanton size distribution. For small- 
sized instantons, such a distribution can be calculated 
by integrating over small gaussian quantum fluctuations 
around the instanton solution. At two loops, 't Hooft 
originally found: 



d(p) 



0.466e- 1679 ^ 



(JV c -l)!(JV t 

Sir 2 
exp 



9 2 (p)J 9 



2)! 
1 



9 2 (P) 



2N C 



(3.4) 



As the size of the instanton increases, the overlapping 
with the field of other pseudo-particles becomes not 
negligible. This generates interaction between pseudo- 
particles which dynamically suppresses the instanton 
weight d{p) for large p. From variational calculations 
[20], phenomenological estimates [lj| and lattice simu- 
lations [21j one finds that the d(p) is peaked around a 
mean value p = 0.3 -7- 0.5 fm. 

In (|3.2|l we have introduced a complex chemical po- 
tential (fi ± i6), associated with the fluctuations of the 
number of instantons and anti-instantons in the ensem- 
ble. Equivalently, one can rewrite the partition function 
in i|3.2|) as a sum over the total topological charge and 
total number of pseudo-particles: 



Z{P; 



From this equation it follows that iO can be interpreted 
as an imaginary chemical potential associated to the fluc- 
tuations of the topological charge. 

The relative width of the fluctuations of N and Q 
around their mean values is proportional to the inverse 
of the size of the system. Thus, in the thermodynamic 
limit, such fluctuations have no consequences on quan- 
tities that are intensive in N and Q. Hence observables 
such as the nucleon mass can be reliably calculated as- 
suming N = (N) and Q — (Q) — 0. On the contrary 
fluctuations can never be neglected when computing av- 
erages of operators which are extensive in N or Q, like 
the case of the correlation function l|3.1|l . 

Since in the ILM topology is clustered around instan- 
tons and anti-instantons, a generic matrix element of the 
type 
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can be written as 12 1 
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where V(Q) denotes the relative occurrence of configu- 
rations with topological charge Q and (0)q is the ex- 
pectation value of the operator O at a fixed topological 
charge Q. Using the fact that \/Q 2 /N — > for large 
TV, Diakonov, Polyakov and Weiss obtained the following 
factorized result: 
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where (Q 2 ) is the topological susceptibility, which in full 
QCD and for small quark masses can be expressed in 
terms of the quark condensate, viz. 
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Fig. 2: Absolute value of the correlation function E T (0, 0) 
as a function of the topological charge Q. Notice that the 
correlator l|2.1l-S|l . which relates directly to the EDM, receives 
contributions only from topologically non-trivial sectors, as 
expected. These results were obtained in the quenched ap- 
proximation, by averaging over the configurations of an en- 
semble of 130 pseudo-particles in a volume V = 3.1 3 x 4 = 120 
fm 4 and quark masses of m q — 0.65 Aqcd- 



This relation holds as well in the instanton vacuum, if the 
fermionic determinant is taken into account |23 , whereas 
in the quenched approximation one has 



(Q 2 ) = N+ + N_ = N 
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Using (|3~lTJ|) and (|3~Tl"|) the neutron EDM in the un- 
quenched case can be expressed as a function of a ratio 
of correlation functions and of the 0-angle parameter: 
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where the quark condensate, the coupling and the mass 
can be computed in the IILM from the appropriate cor- 
relation functions. 

Similarly, in the quenched approximation the electric 
dipole moment reads 
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Notice that the above Eqs. (|3.13|l and l|3.14|l are well 
defined in the thermodynamic limit, because the deriva- 
tive of E T is proportional to the inverse volume and the 
number of pseudo-particles N grows linearly with the 
volume. 



4. CHIRAL BEHAVIOR OF THE NEUTRON 
EDM WITHIN THE IILM 

In this Section we derive the qualitative behavior of 
the expressions (|3.13|) and Q3.14J I in the chiral limit by 
studying the two quantities appearing in the r.h.s. of 
Eq. (|3.10(1 . We will not consider corrections arising from 
chiral loops. 

Let us start with the unquenched case, where the topo- 
logical susceptibility (Q 2 ) vanishes linearly in the quark 
mass 
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Let us now look at the second factor in the r.h.s. of 
Eq. (|3.10() . Writing out the discrete derivative explicitly 
one has 
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where we have used the fact that the average of our CP 
violating operator vanishes in trivial topological sectors 
(Q = 0). Writing the average (0)q explicitly in terms of 
functional integrals one gets 
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where TrfS'S'...] denotes a trace of propagators arising 
from the explicit integration over the fermionic fields 
(Wick contractions). 

In order to study the chiral behavior of the trace we 
recall that a quark propagator in the given background 
can be written as: 
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In a topologically non-trivial sector of the instanton 
vacuum, the Dirac operator has Q exact zero-modes as- 
sociated to the extra instantons (or anti-instantons) in 
the ensemble: 

$4=0, j = l,...,Q. (4.6) 

The quark propagator can then be written as: 

s(x y) = y^ ^JM^oli^ | y^^i^liy) 



-im f 



3 = 1 J A^O 

^>top\X, y) + Onon~top{X, 2/) 



A — irrif 



(4.7) 



In this expression S top (x,y) is the topological contribu- 
tion to the propagator arising from the exact zero-modes. 
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Fig. 3: ILM results obtained at different values of the quark 
masses with quenched (circles) and unquenched (squares) sim- 
ulations. The behavior of the unquenched and quenched cal- 
culations for small quark masses is consistent with Eqs. (14.91 
and (14.101 . respectively. 



It is important to recall that the Pauli principle implies 
that at most Nf quarks can propagate in the exact zero- 
mode wave function of an instanton. Thus one has 
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Collecting the results for chiral behavior of all the 
terms relevant for the EDM, one gets 
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Hence we conclude that, when the fermionic determi- 
nant is included, the neutron EDM vanishes linearly in 
the quark mass, as it is also expected from other es- 
timates based on chiral perturbation theory (see, e.g., 
Ref. 0). 

When the quenched approximation is adopted, the 
topological susceptibility (Q 2 ) becomes independent of 
the quark mass, while neglecting the fermionic determi- 
nant removes a factor m f from the numerator of the 
r.h.s. of Eq. (|4.9(l . Therefore one expects a divergent be- 
havior in the chiral limit of the form 



D M 



1 



m N f 



(4.10) 



We stress the fact the mass dependencies given in 
Eqs. (|4.9(l and l|4.10|l do not depend on the particular 
values of the model parameters which define the ILM. 
These results rely only on the working assumption that 
the quantum mixing of the QCD vacuum can be de- 
scribed in terms of isolated tunneling events (instantons). 
In such a semi-classical limit, the Index Theorem is re- 
alized in a very specific way, with Nf exact zero-modes 
associated to each unit of topological charge. 



V = 56 fm 4 
m = 1.00 Aqcd 


|W>| 1/3 [GeV] 


M N [GeV] 


Z [GeV 6 ] 


quenched 


0.186 


1.42 


0.018 


unquenched 


0.186 


1.41 


0.017 




V = 90 fm 4 
m = 0.75 Aqcd 


|W>| 1/3 [GeV] 


M N [GeV] 


Z [GeV 6 ] 


quenched 


0.199 


1.32 


0.018 


unquenched 


0.187 


1.30 


0.016 




V = 120 fm 4 
m = 0.65 Aqcd 


|(V^>| 1/3 [GeV] 


M N [GeV] 


Z [GeV 6 ] 


quenched 


0.204 


1.27 


0.018 


unquenched 


0.199 


1.10 


0.012 



Table 1: Results of quenched and unquenched simulations 
for the quark condensate (ip4>) , the neutron mass M and 
the constant Z. The parameter Aqcd is taken at the value 
Aqcd = 220 MeV. 



V = 56 fm 4 
m = 1.00 Aqcd 


(Q 2 ) 


\dT, T /dQ\ [GeV 8 ] 


D/6 [e • fm] 


quenched 


56 


0.015 


0.35 ±0.06 


unquenched 


22 


0.025 


0.32 ±0.04 




V = 90 fm 4 
m = 0.75 Aqcd 


(Q 2 ) 


\dT, T /dQ\ [GeV 8 ] 


D/6 [e • fm] 


quenched 


90 


0.043 


0.79 ±0.07 


unquenched 


32 


0.039 


0.27 ±0.07 




V = 120 fm 4 
m = 0.65 Aqcd 


(Q 2 ) 


\dZ T /dQ\ [GeV 8 ] 


D/6 [e • fm] 


quenched 


120 


0.064 


1.04 ±0.06 


unquenched 


40 


0.046 


0.17 ±0.04 



Table 2: Results of quenched and unquenched simulations 
for the topological susceptibility (Q 2 ), the quantity dT, T /dQ 
and the neutron EDM divided by the angle 6. The statistical 
error in the last column is dominated by the uncertainty on 
the quantity dT, T /dQ. 



5. NUMERICAL RESULTS 

We have performed quenched and unquenched numer- 
ical simulations in the IILM, assuming two active de- 
generate flavors. We have used three values of the quark 
mass m, namely m = 0.65, 0.75, 1.0, in units of the Pauli- 
Villars scale parameter Aqcd = 220 MeV [l]|. Corre- 
spondingly, we have used three periodic boxes of volume 
given by V = (3.1 3 x 4), (2.8 3 x 4) and (2.4 3 x 4) fm 4 . 
The total density of pseudo-particles in these ensembles 
has been chosen to satisfy (N + ± N-)/V = 1 fm -4 . 

We have computed the quantities which are needed 
to determine D/6 through Eqs. l{3~T3j) and Qttty . The 
neutron mass M and the coupling Z of the interpolating 
operator have been obtained by extracting the slope and 
intercept of the logarithm of the two-point function (see 
the example in Fig. QJ. The derivative of (|2.13|l with 
respect to Q has been computed by varying by a small 
amount the number of instantons relative to that of anti- 
instantons in the box, while keeping the total number of 
pseudo-particles fixed (see the example in Fig. [5J ■ 

The results of our calculations are summarized in Ta- 
bles H21 From Table Q] we observe that, in case of the 
quark condensate, the neutron mass and the constant Z 7 
quenched and unquenched simulations give almost the 
same results. This is as expected, since these quantities 
are not directly related the topological structure of the 
QCD vacuum. On the contrary, the quantities reported 
in Table[5]show a much larger sensitivity to the effects of 
the inclusion of the fcrmionic determinant. In particular 
the topological susceptibility (Q 2 ) [see Eqs. ( I3.11|) and 
lj3~l"!>|) ] and the EDM D/6 [see Eqs. (J3~T3|) and IjiTHjl ] 
can differ by a factor of ~ 3 ± 4 between the quenched 
and unquenched simulations with our choice of the quark 
mass values. 

The quenched and unquenched results obtained for the 
EDM are also shown in Fig. [3] in order to better appreci- 
ate the quark mass dependence. The unquenched results 



appear to be consistent with a linear dependence on m 
with a zero intercept, as expected from QCD and from 
the arguments leading to Eq. Ij4.9|) . At the same time the 
quenched results exhibit a sharp increase at low quark 
masses, consistent with a divergence of the form given 
by Eq. (|4.10() with Nf — 2. Clearly, a quenched calcula- 
tion of the neutron EDM becomes completely unreliable 



for 



< A 



QCD- 



In the quenched approximation, the divergence appear- 
ing in the chiral limit makes it impossible to perform the 
extrapolation toward the physical quark mass value. On 
the other hand, such an extrapolation is possible in the 
case of unquenched results. Assuming a linear mass de- 
pendence, the IILM prediction for the neutron EDM cor- 
responding to a light quark mass between 4 and 10 MeV 
is: 



|d„| = (6 ±14) x 10~ 3 6 (c • fm) 



6. CONCLUSIONS 



(5.1) 



We have used the Instanton Liquid Model to study the 
role played by the fcrmionic determinant in the evalua- 
tion of the neutron EDM, which is an observable sensitive 
to the topological structure of the vacuum. We have ana- 
lyzed the chiral behavior of such a quantity (up to chiral 
logs) both in the quenched and unquenched cases. We 
have found that, when the fcrmionic determinant is in- 
cluded, the neutron EDM is expected to vanish linearly 
with the quark mass, whereas in the quenched approxi- 
mation it should diverge as 1/m f in the chiral limit. 

We have performed several model simulations and 
found that quenched and unquenched calculations give 
comparable results for the neutron EDM at large quark 
masses (~ Aqcd), whereas they strongly differ at lower 
quark masses. At the lowest value of the quark mass used 
in our simulations (m ~ 130 MeV) the quenched result 



8 



is a factor of ~ 4 larger than the unquenched one. 

We have obtained the ILM prediction for the neutron 
EDM by extrapolating the unquenched IILM result to 
the physical value of the quark mass. The ILM result is 
roughly a factor 2—4 larger than existing model estimates 

m. 

Our main conclusion is that quenched and unquenched 
lattice QCD simulations of the neutron EDM as well as 
of other observables governed by topology might show 
up similar important differences in the quark mass de- 
pendence, near the chiral limit. In particular, our semi- 
classical analysis suggests that a quenched lattice calcula- 
tion of the neutron EDM could be affected by a topology- 



driven divergence, which would make it impossible to per- 
form the extrapolation to the physical value of the quark 



We insist on the fact that the qualitative predictions 
in Eq. (|4.9(1 and (|4. 10|> do not depend on the particular 
values of the model parameters which define the ILM. 
They rely only on a semi-classical description of the quan- 
tum mixing in the 6- vacuum, in terms of isolated tunnel- 
ing events. Hence, the observation of a divergence in a 
quenched lattice calculation of the neutron EDM, would 
represent a clean, parameter-free signature of instanton- 
induced dynamics in QCD. 
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